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We derive and solve analytically the non-Markovian master equation for harmonic quantum Brow- 
nian motion proving that, for weak system-reservoir couplings and high temperatures, it can be re- 
cast in the form of the master equation for a harmonic oscillator interacting with a squeezed thermal 
bath. This equivalence guarantees preservation of positivity of the density operator during the time 
evolution and allows one to establish a connection between the dynamics of Schrodinger cat states 
in squeezed environments and environment-induced decoherence in quantum Brownian motion. 
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I. INTRODUCTION 

The increasing ability in coherent control and manipu- 
lation of the state of quantum systems has paved the way 
to experiments that are able to monitor the transition 
from quantum superpositions, such as the Schrodinger 
cat states, to classical statistical mixtures [J The 
emergence of the classical world from the quantum world, 
due to decoherence induced by the environment, has been 
extensively investigated in the last few decades both in 
connection to the fundamental issues of quantum theory 
and in relation to the emerging quantum technologies. 
The fragile nature of quantum superpositions and entan- 
gled states, which are exploited, e.g., in quantum commu- 
nication, quantum computation, and quantum metrol- 
ogy, makes these potentially very powerful techniques 
also very delicate j^]. 

In order to understand the quantum to classical tran- 
sition in terms of decoherence induced by the environ- 
ment a paradigmatic model of the theory of open quan- 
tum systems, namely the quantum Brownian motion 
(QBM) model, has been extensively used in the litera- 
ture I, S m B i, M HH, M M, 0, d, M M M- 

This model can be solved exactly using the path inte- 
gral approach. This numerical method involves, however, 
complicated integrals and therefore simpler analytical ap- 
proaches are desirable. In the weak-coupling limit, a sec- 
ular master equation has been derived and investigated 
extensively in Refs. [ul, [H, [H [ll] ■ The secular ap- 
proximation simplifies greatly the study of the dynamics 
and allows one to understand the origin of microscopic 
physical processes, such as the virtual exchanges of en- 
ergy between the system and the environment, character- 
izing the non-Markovian dynamics. This approximation 
can be safely used when dealing with a class of observ- 
ables that, in the weak coupling limit, do not depend 
on the non-secular terms, such as the mean energy of 
the system [ll], [l^. The study of the non-Markovian 
short time behavior of the density matrix of the quan- 
tum Brownian particle, however, requires in general the 
use of the complete QBM master equation. 

It is well known that the weak coupling limit of the 
QBM master equation is operatorially not in the Lind- 
blad form [s^]. This fact complicates the analysis of the 



system dynamics and requires one to resort to numerical 
approaches. Moreover, the positivity of the density ma- 
trix during the time evolution can be violated, leading to 
unphysical results [l9|. However, by adding some terms 
that are negligible in the weak coupling limit, the Marko- 
vian master equation can be rewritten in the Lindblad 
form (see, e.g., Ref. Q). The main result of this paper 
is to establish a connection between the weak coupling 
master equation for QBM and the master equation for a 
harmonic oscillator interacting with a squeezed thermal 
reservoir [13, HH, [H, HI, HI. The latter one has been 
extensively studied in the quantum optical context, and 
the conditions under which the positivity of the density 
matrix is preserved during the time evolution are known. 
Using the connection between these two master equations 
we obtain an analytic solution to the quantum Brown- 
ian motion dynamics for a weak coupling but in pres- 
ence of a structured reservoir, i.e., without performing 
the Markovian approximation. Our analytic solution al- 
lows one to give a simple and physically clear description 
of environment-induced decoherence of a Schrodinger cat 
state for the QBM model, and provides results in agree- 
ment with Ref. 16], in the Markovian limit. 



Very recently the quantum Brownian motion model 
has been used to study the entanglement between two 
harmonic oscillators in a common bath 25, 26, 27, 28]. 
In this context a very rich scenario showing the appear- 
ance of entanglement sudden death and revivals has been 
brought to light. We believe that our results may shed 
light on the microscopic origin of these phenomena. 



The outline of the paper is the following. In Sec. II we 
recall the master equation for quantum Brownian mo- 
tion and its properties. In Sec. Ill we derive the weak 
coupling non-Markovian master equation and we estab- 
lish the connection with the master equation describing 
a harmonic oscillator in a squeezed thermal bath. Sec- 
tion IV is devoted to the solution of the master equation 
and to the study of the non-Markovian dynamics for an 
initial Schrodinger cat state. Finally, Sec. V contains the 
conclusions. 
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II. THE QBM MASTER EQUATION 

We consider a harmonic oscillator linearly coupled to 
an engineered reservoir modelled as an infinite chain of 
non-interacting oscillators 0, H, 0, 0, H, H, The mi- 
croscopic Hamiltonian of the total system, in units of h, 
has the form 



uJo[a^a + -] , 



1 



H = Ho + He + Hi 
Ho 

He 

Hint 



g^i^n {a'' + a) + b„ 



(1) 
(2) 

(3) 

(4) 



where Hq, He, and Hi^t are the system, environment 
and interaction Hamiltonians, respectively, a (a^) and 
bn are the annihilation (creation) operators of the 
system and of the reservoir quantum oscillators, respec- 
tively, luq is the frequency of the system oscillator, a;„ 
are the frequencies of the reservoir oscillators, g is the 
coupling constant, and the quantities k„ describe how 
strongly the reservoir oscillators are coupled to the sys- 
tem. In the continuum limit one introduces the spec- 
tral distribution J{uj) related to k„ via the equation 
J{uj) — J2n i^nS{uj — cj„)/(2m„a;„), with to„ as the mass 
for each environmental oscillator 5]. 

Following the standard derivation of the master equa- 
tion for the reduced system (see, e.g., Ref. @) one can 
demonstrate that the exact master equation, in the in 
teraction picture, takes the form 



ier equation 



dpjt) 
dt 



+ m[xAp.pm + \r{t)[x\p{t)] 



(5) 



where p[t) is the reduced density matrix, X = 
(a + a^) 1^/2 and P = i (a'f - a) /\/2. 

The master equation given by Eq. ([5]), being exact, 
describes also the non-Markovian short time system- 
reservoir correlations due to the finite correlation time 
of the reservoir. In contrast to other non-Markovian dy- 
namical systems, this master equation is local in time, 
i.e., it does not contain memory integrals. All the 
non-Markovian character of the system is contained in 
the time-dependent coefficients appearing in the master 
equation (for the analytic expression of the coefficients 
see, e.g., Ref. [13j). These coefficients depend uniquely 
on the form of the reservoir spectral density. The co- 
efficient r{t) describes a time-dependent frequency shift, 
7(t) is the damping coefficient, A(i) and n(t) are the nor- 
mal and the anomalous diffusion coefficients, respectively 

The dynamics of the system has been extensively stud- 
ied numerically using the path integral approach (see 



Refs. [1, for a review). In particular, this model 
has been used to demonstrate the action of environment- 
induced decoherence for initial Schrodinger cat states 
such as, e.g. 



I*) 



where |a) is a coherent state. 



AA-i =2[l + exp(-2|ap)] , 



(6) 



(7) 



and we take a G M for simplicity. It has been proven 
that the decoherence induced by the environment acts in 
a much faster time scale than the thermalization process, 
in particular the decoherence time Td is inversely pro- 
portional to the separation between the two components 
of the superposition, i.e. 2|ap. For this reason, when 
the environmental spectrum has a structure, the reser- 
voir correlation time, characterizing the duration of the 
non-Markovian dynamics, can become comparable to the 
decoherence time scale. In this case a non-Markovian de- 
scription of environment-induced decoherence is impor- 
tant even in the weak coupling limit. 



III. CONNECTION WITH THE MASTER 
EQUATION FOR A SQUEEZED THERMAL 
BATH 

For the sake of definiteness we focus on an Ohmic reser- 
voir described by a spectral distribution of the form [3| 



J{u) 



2w 



(8) 



with as the cutoff frequency. The extension of our 
results to other forms of spectral densities is straight- 
forward. In the limit of high temperatures and for suffi- 
ciently weak system-reservoir couplings, both the anoma- 
lous diffusion term n(t) and the frequency shift term r{t) 
are negligible [111, ^3] and the QBM master equation @ 
can be recast in the form 



dp{t) A(t)+7(0 ^, t^ , m-l{t) r^ 1 
^ ^ 2 ) + 2 



+ 2 



A(i)_ 



2iuJot 



D[a\ (9) 



where the superoperators L[0\ and D[0] are given by 

L[0] = OO^p + pOO^ -20^pO, 

D[0] = 0'^p + pO^ -20pO. (10) 

Equation ([9]) is not in the Lindblad form, but has the 
same operatorial structure of the master equation for 
a harmonic oscillator in a thermal squeezed bath. In- 
deed, having in mind, e.g., Eq. (6) of Ref. [l^ we can 
see that the two master equations coincide provided that 
we take A(t) + -fit) = -f{N + 1), A{t) - 7(t) = 77V 
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and A(t)e2^'^°* = -7M. The condition for positiv- 
ity of the density operator, for a squeezed reservoir, 
reads as follows |Mp < N{N + 1). It is straightfor- 
ward to verify that Eq. Q violates this condition since 
|M|2 = A{t)^ > A{tf - 7(t)2 = N{N + 1). 

In order to cure this problem we take a closer look at 
the second order expansion of the diffusion and dissipa- 
tion coefficients appearing in Eq. ^ 10, 13] 

poo 

A{t)= g'^ ldujdtiJ((J)[2N{uj) l]cos(tjti)cos(woii), 

(11) 

nt poo 

j{t) = dwdii J(a;) sin(ti;ti) sin(tJo^i), (12) 

Jo Jo 

where N{uj) — {^c^/^bT _ 1)^^ is the average number of 
reservoir thermal photons, fcs is the Boltzmann constant, 
and T is the reservoir temperature. We note that, for 
high T, i.e., N{io) > 1, we have A{t) > 7(t). 

Having this in mind we can modify the weak coupling 
master equation ([9]) by adding some terms, proportional 
to 7(t), which are neghgible with respect to the terms 
proportional to A{t) in the high T and weak couphng 
limits. The resulting master equation takes the form 



dp{t) 
dt 



m+lit) r> t^ , ^(^)-7(^) ,r 1 



D\a' 



(13) 



In this case A{t) + 7(t) = -f{N + 1), A{t) - -/{t) = -/N 



and [A{t) - j{t)]e 



2iujQt — 



'^M, and the condition for 



positivity of the density operator |Afp < N{N + 1) is 
satisfied as an equality. This corresponds to the case of 
a maximally squeezed reservoir (|M| = N). 

The analytical form of the coefficients A{t) and ^{t) 
can be easily calculated explicitly for an Ohmic environ- 
ment, in the high T limit, inserting Eq. ([8]) into Eqs. 



A{t) = 2g^kBT- 



1 + 

- (1/r) sin(c^oi)]} 



{1- 



[cos(Li;o^) 



(14) 



lit) 



2 2 

g ujor^ 



1 + 



1-6"""=* cos(woO- sm{uJot) 



(15) 



with r ~ LUc/ujQ. 

From these equations we see that A{t) and j{t) start 
from an initial zero value and quickly approach their con- 
stant Markovian value. Indeed for tr = 1/wc, with 
Tfi as the reservoir correlation time, the time-dependent 
coefficients A(t) + ^{t) and A{t) — ^{t) become 



m+i{t) 

A{t) - 7(i) 



t[n{loo) + i] = ri 

VN{LOa) EE 7^i 



M 



(16) 
(17) 



respectively, with N^loq) ~ ksT/LOQ and 



29' 



+ 1 



LOq. 



(18) 



In the next section we will present the solution of this 
master equation ()13p in terms of the Wigner function 
and discuss its properties for an initial state of the form 
of Eq. Moreover, we will check the validity of the ap- 
proximations under which this master equation holds by 
comparing the solution we have derived with the solution 
of the exact Hu-Paz-Zhang master equation ([5]) . 

By looking at Eq. p3)) one realizes immediately that, 
for times loqI ^ 1 the last two terms average out to zero 
and the master equation reduces to the secular master 
equation used, e.g., in Refs. [H, [13, [H, HO] 



dpjt) ^ A{t)+^{t) 
dt 2 



[2ap{t)a' — a) ap{t) — p{t)a) a\ 



\2a) p{t)a — aa) p{t) — p{t)aa)"\ 



(19) 



It is known that in the weak coupling limit there exists 
a class of observables, e.g. n — a^a, whose dynamics is 
not affected by the counter-rotating terms, i.e. by those 
terms neglected in the secular approximation |Tl[. 
In general, however, the counter-rotating terms do con- 
tribute to the dynamics of the reduced density operator 
of the system, in particular in the short non-Markovian 
time scale we are interested in. 

In the following we will focus on two specific physical 
regimes characterized by opposite values of the parame- 
ter r, namely the case in which r — Wc/wq ^ 1 and the 
case r ^ 1. We will refer to these cases as the resonant 
and the off-resonant case, respectively, since for r ^ 1 
the frequency cjq of the oscillator overlaps with the spec- 
trum of the reservoir while, for r <C 1, wo is off-resonant 
with the spectrum J{i-u). From previous studies on the 
dynamics of the heating function {n{t)) we know that 
the system time evolution differs notably in these two 
regimes As we will see in Sec. IV, also the dynam- 
ics of an initial Schrodinger cat state crucially depends 
on the value of the resonance parameter r. 

It is worth emphasizing that, in the resonant case r 3> 
1, for ujct < 1, we cannot use the secular master equation 
since ojot <C oJct < 1. Therefore we should use Eq. (|13p . 
On the contrary, in the off-resonant case r ^ 1, we can 
focus on the dynamics for times 1/ujq <^ t < l/uc, since 
this is consistent with the assumption r — lUc/ujq <^ 1, 
and use the simpler secular master equation (fTH)) . 



IV. THE DYNAMICS 

We now look at the short-time dynamics of a 
Schrodinger cat state of the form given in Eq. ©. This 
state is also known as even coherent state due to the fact 
that only the even components of the number probability 
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FIG. 1: (Color on line) Wigner function W{l3) and number 
probability distribution P„ for the state given by Eq. ^ with 
a = 2. 



distribution are nonzero. The oscillations in the number 
state probability are a strong sign of the nonclassicality of 
this state. This and other nonclassical properties of the 
even coherent state, such as the negativity of the corre- 
sponding Wigner function, have been extensively stud- 
ied in the literature (see, e.g., Ref. [20| and references 
therein). In Fig. 1 we show the number probability dis- 
tribution and the Wigner function for the even coherent 
state. This state has been realized in the trapped ion 
context and the transition from a quantum superposi- 
tion to a classical statistical mixture has been observed 
experimentally 

The decoherence and dissipation due to the interaction 
with both thermal reservoirs and squeezed reservoirs has 
been studied, in the Markovian limit, in Refs. [l^l and 
pll . [2^ . [23I I , respectively. On the other hand, the exact 
dynamics of the Hu-Paz-Zhang master equation ([S]) has 
been numerically investigated in Refs. [3, [IB]- In the 
following we will look for a simple analytic solution valid 
in the short-time non-Markovian regime in the two limits 
r <C 1 and r 3> 1. We will verify the validity of the 
approximations made to derive such a solution comparing 
it with the Hu-Paz-Zhang result 



In order to describe the transition from the initial even 
coherent state to the corresponding classical statistical 
mixture induced by the interaction with the environment 
it is convenient to look at the dynamics of the Wigner 
funtion. In order to investigate the non-Markovian dy- 
namics we need to find out the time evolution of the 
Wigner function for times t < I/lUc- 



1. The off-resonant case r ^ 1 

In the off-resonant case r <C 1 the solution of Eq. (fTQ]) 
in terms of the quantum characteristic function (QCF) 
x(^) is presented in Ref. [ll|. From this solution, and re- 
membering that the Wigner function is the Fourier trans- 
form of the QCF, 

1 r'^ 

W{p) = — / £ix{0 exp {PC - P*0 , (20) 
^ J —00 

one obtains straightforwardly the following expression for 
the dynamics of an initial even coherent state. 



WiP,t) =W^+"\p,t) +W^-''Hl3,t) + WiiP,t), (21) 
with 

^f ( 01 



TT \N<yt) + 1/2] 



X exp 



exp 



,-r(t)/2. 



iV(i) -f 1/2 



N(t) -I- 1/2 



(22) 



2M 



TT [N{t) + 1/2] 



exp 



N{t) + 1/2 



X exp 



-2a' 1 



2N{t) + 1 



2e-r(t)/2 
N{t) + 1/2 



with 



N{t) 

m 



dt'A{t'), 



dt'i{t'), 



(23) 

(24) 
(25) 



where the coefficients A(t) and 7(t) are given by 
Eqs. (|14p - p3|) . As known from the Markovian theory, 
the interaction with the environment causes the disap- 
pearance of the interference peak and therefore the tran- 
sition from quantum superposition to classical statistical 
mixture. A useful quantity to monitor this transition is 
the fringe visibility function 



F{a,t) = exp(-A,„f) 

^ 1 ^/(/?.^)lpeak 



(26) 



peak 



r/2 ■ 



where we indicate with i) I peak 



peak 



the value of the Wigner function 



at /3 = (0,0) and (3 — (±Q!,0), respectively. Inserting 
Eqs. (Illll-dlSI) into Eq. ^ we obtain, for the r < 1 
case, 



F{a, t) — exp 



-20" ( 1 



-r(t) 



2N{t) -f 1 



(27) 
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This equation tells us that, as for the Markovian the- 
ory, the interference term disappears faster and faster, 
the larger is the separation between the two components 
of the superposition, measured by 2q!, and the higher is 
the temperature of the environment [see Eqs. (|14p - ((24|) ]. 
In Fig. 2 we plot the time evolution of the fringe visibil- 
ity factor comparing the off-resonant (dashed line) and 
the on-resonant (solid line) cases. 



Z. The resonant case r S> 1 




rt 



We now consider the dynamics in the more complicated 
resonant case. We begin by noting that, for ujct < 1, we 
can approximate 6^*"^°* ~ 1 since r = oJc/^o ^1- As a 
consequence we have M S R and M = —N. 

The solution of the master equation p3|) can be ob- 
tained following the same steps as in solving the mas- 
ter equation for a harmonic oscillator in a thermal 
squeezed bath, in the Markovian case. We write down 
the Fokker-Planck equation corresponding to the master 
equation for the Wigner function. Since the initial 
state is a linear combination of Gaussian terms, the form 
of the master equation ensures that each Gaussian term 
evolves independently. Therefore the evolved state will 
also be a linear combination of Gaussian terms [13] . Hav- 
ing this in mind it is straightforward to derive the Wigner 
function dynamics as follows 



W±"(/3,t) = 



Tr{N{t) + l/A}^/^ 



X exp 



X exp 



2N{t) + l/2 



-r(t)/2^ 



1/2 



(28) 



Wiif3,t) 



2N 



X exp 



X exp 



X cos 



2N{t) + l/2 1/2 
-2a' ( 1 - 



AN{t) + 1 



2e-r(t)/2 
2N{t) + 1/2 



a[3i 



(29) 



with N{t) and V{t) given by Eqs. ([24]) -([25]). Comparing 
the dynamics to the case r <C 1, where the secular ap- 
proximation holds, we see the following differences. First 
of all, the variances of the two Gaussians M^^"(/3, t), cor- 
respondent to the two components of the superposition, 
do not follow the same dynamics anymore. This asym- 
metry in the time evolution of the variances is typical 
of the behavior of harmonic oscillators in squeezed en- 
vironments. But more interesting for the study of the 
quantum-classical transition is the behavior of the fringe 



FIG. 2: (Color on line) Time evolution of the fringe visibility 
factor for the state given by Eq. (O with a = A, kBT/cuo = 
100, r = 10 (solid line) and r = 0.1 (dashed line). 



visibility 

i^(a, t) — exp 



-2a' 1 - 



.-r(t) 



AN{t) -f 1 



(30) 



Comparing this equation with Eq. ((27)l obtained for 
r ^ 1 we notice that, with regard to the decoherence 
process, in the resonant case r 3> 1 it is as if the system 
would interact with a thermal reservoir with an effective 
temperature that is the double of the real temperature. 
This is simply due to the non-negligible role played by 
the counter-rotating terms present in the microscopic in- 
teraction Hamiltonian of Eq. ^ . In Fig. 2 we plot the 
time evolution of the fringe visibility factor for an initial 
Schrodinger cat state with a = 4, in the cases r = 10 
(solid line) and r = 0.1 (dashed line). 

To conclude this section we compare the dynamics of 
the fringe visibility factor obtained without performing 
the secular approximation, i.e. by solving Eq. (jlSp . with 
the result presented in Ref. [11] for the exact Hu-Paz- 
Zhang master equation ([5]). Combining Eq. ((26|) and 
Eq. ([201), and using Eq. ([Mj) we obtain 



2a' 



^4/Jdi'A(i')4 



1 



(31) 



where we have put e"'"^*-' ~ 1 since we are far from 
thermalization, i.e., t ^ tth- We recall that after a 
time t ~ Tfl the diffusion coefficient attains its constant 
Markovian value Am and therefore /p dt'A{t') ~ Am^- 
Having this in mind we can compare Eq. pip with the 
high temperature and low damping approximation of 
Aint given by Eq. (42) of Ref. 16]. One sees immedi- 
ately that, mutatis mutandi, our Eq. pip coincides with 
Eq. (42) of Ref. [l^ for t ^ th and also it provides 
a straightforward generalization of this equation for the 
case t <tr. 

Another important aspect, discussed in Ref. [3] and 
predicted by the exact Hu-Paz-Zhang model, is the differ- 
ent decoherent dynamics of superpositions of wave pack- 
ets separated in positions, such as those we considered 
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in our example, and of superposition of wave packets 
separated in momentum. The latter ones experience a 
much slower decoherence compared to the former ones. 
In Ref. this fact was attributed to the microscopic 
interaction Hamiltonian coupling the position operator 
X of the system oscillator to the bath. This prediction 
is confirmed by our model. For times t <C tthorm, with 
^thorm the thermalization time, indeed, we can neglect 
the small difference between the rates A(t) + ^{t) and 
A{t) — 7(t), corresponding to the upward and downward 
transition, i.e., to the absorption and emission of a quan- 
tum of energy by the system oscillator, respectively. In 
this case the master eqution becomes 



dpit) Ait) ^ Ait) 



Ait)[X, [X,p]], 



(32) 



describing a continuous measurement of the position of 
the system harmonic oscillator. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied the non-Markovian short 
time dynamics of the harmonic QBM model. In the high 
temperature case we have obtained a weak coupling mas- 
ter equation preserving positivity of the density matrix 
and able to describe situations where the spectrum of the 



environment is not flat. This master equation has the 
same form as the master equation describing a harmonic 
oscillator interacting with a squeezed thermal reservoir 
and preserves positivity of the density matrix during the 
time evolution, thus solving the problem discussed, e.g., 
in Ref. ;19]. 

Using this analogy, we have obtained a solution in 
terms of the Wigner function and investigated the dy- 
namics in two opposite physical regimes characterized by 
the parameters r ^ 1 (resonant case) and r ^ 1 (off 
resonant case). A comparison between the two solutions 
reveals that the counter-rotating terms present in the mi- 
croscopic Hamiltonian model have a non-negligible effect 
in the resonant case. On the contrary, in the off-resonant 
case they can be neglected. 

The availability of simple analytic expressions describ- 
ing the short time dynamics of the system oscillator has 
allowed us to shed light on the non-Markovian dynamics 
of environment induced decoherence and, in particular, 
on its sensitivity to the form of the initial superposition. 

Because of the fragility of such states it is interesting 
to see whether it is possible to modify the decoherence 
induced by the environment, e.g., by means of the quan- 
tum Zeno effect [sO] • We plan to investigate this issue in 
the future. 
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